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Abstract — In ultra-wideband (UWB) communication systems 
with impulse radio (IR) modulation, the bandwidth is usually 
1GHz or more. To process the received signal digitally, high sam- 
pling rate analog-digital-converters (ADC) are required. Due to 
the high complexity and large power consumption, monobit ADC 
is appropriate. The optimal monobit receiver has been derived. 
But it is not efficient to combat intersymbol interference (ISI). 
Decision feedback equalization (DFE) is an effect way dealing 
with ISI. In this paper, we proposed a algorithm that combines 
Viterbi decoding and DFE together for monobit receivers. In this 
way, we suppress the impact of ISI effectively, thus improving 
the bit error rate (BER) performance. By state expansion, we 
achieve better performance. The simulation results show that 
the algorithm has about ldB SNR gain compared to separate 
demodulation and decoding method and ldB loss compared to 
the BER performance in the channel without ISI. Compare to 
the full resolution detection in fading channel without ISI, it has 
3dB SNR loss after state expansion. 

Index Terms — Monobit, decision-feedback equalization, joint 
decoding, ultra-wideband 

I. Introduction 

Impulse radio ultra-wideband (IR-UWB) systems utilize the 
short pulse with a low duty cycle to carry information [1 1. The 
short duration of the pulse enables high data rate but occupies 
large bandwidth, usually 1 GHz or more. In communication 
systems, digital approaches can provide attractive flexibility 
in receiver signal processing. In order to process the received 
signal digitally, the analog-digital-converter (ADC) is required. 
Due to the ADC resolution and the limitation of its power 
consumption, it is appropriate to use low-resolution ADC. 
Monobit ADC has a simple structure and can be realized by a 
fast comparator. Tens of Giga samples per second (Gsps) could 
be reached. Monobit ADCs have low power consumption and 
low cost. 

Monobit digital receivers for IR-UWB have been developed, 
see e.g. 0. Optimal monobit receiver under Nyquist sampling 
rate has been proposed in Q. The optimal receiver turns out 
to be a linear combiner. Suboptimal monobit receiver was also 
proposed. 

In indoor channel, due to the multipath propagation, the 
communication system suffers from the effect of intersymbol 
interference (ISI). The receiver has to detect symbols in 
the ISI. The optimal monobit receiver mentioned above was 
derived under the assumption that maximum channel delay is 



significantly smaller than symbol duration. Thus it is may not 
perform well when the delay spread is large. 

There are plenty of ways to cope with ISI H). Current 
strategies that suppress ISI mainly include the following: the 
equalizer based on minimum mean square error (MMSE) rule, 
which is introduced and analyzed in 0, the Zero Forcing (ZF) 
Equalizer in the Bell System Technical Journal in 1965 [6| and 
the decision-feedback equalizer (DFE) which cancels the ISI 
by feeding back the decided results Q. Among them, DFE is 
relatively a simple approach. At high SNR, this method makes 
reliable decisions. However, at low SNR, error propagation 
is a severe problem. In practical communication systems, 
channel coding is a critical component. Under full resolution 
(FR) sampling, a joint coding and DFE method has been 
proposed in iflOll . This method uses a combination of Viterbi 
soft decisions and delayed decisions to minimize bit error rate 
(BER). Unfortunately, in UWB systems, as a result of the 
short symbol duration, full resolution received waveform is 
difficult to obtain. For monobit receivers, no algorithms have 
been developed to combat heavy ISI as far as we know. 

In this paper, we propose a joint Viterbi decoding and DFE 
algorithm for monobit receivers. This algorithm provides an 
efficient way for monobit receivers dealing with ISI in indoor 
fading channel. The channel has one line-of-sight (LOS) path 
and many non-line-of-sight (NLOS) paths. The NLOS paths 
appear in a relatively long time after the LOS arose. In 
this algorithm, for every arrived state in the code trellis, we 
construct a full-resolution reference waveform to cancel the 
overlapped ISI, then we use derived monobit optimal approach 
to compute the path metric f2| ■ The path having the maximum 
likelihood probability metric remains and is set as the current 
surviving path. We will derive the likelihood probability for 
each state. State expansion is also proposed to improve the 
BER performance. The simulation results show that compared 
to the fading channel without ISI, the receiver suffered 2dB 
SNR loss from ISI. After state expansion, the loss reduces 
to ldB. Compared to FR detection in the fading channel 
without ISI, the receiver has 3dB SNR loss. Compared to 
separate demodulation and decoding, the joint decoding and 
DFE algorithm has ldB SNR gain. 

The rest of the paper is organized as follows: In Section 
II, we describe the system model. Section III describes the 
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Fig. 1. System Blockdiagram 



algorithm we proposed. We give an example and discuss state 
expansion method. We can use this approach get better BER 
performance but the complexity increases. Simulation results 
are provided in Section IV. Section V is the conclusion. 

II. System Model 

The block diagram of joint Viterbi decoding and DFE 
monobit digital receivers is depicted in Fig. Q] The baseband 
received signal is first filtered by an ideal low pass filter (LPF). 
The bandwidth of the filter is B. Then the received signal will 
be sampled at Nyquist sampling rate T — 1/2B and quantized 
to one bit resolution. The digitized signal is processed by a 
digital signal processing (DSP) unit for decoding and symbol 
estimation. 

In this paper, the transmitted information are binary sym- 
bols. We assume every data block has U binary symbols. 
d u € {+1,-1} is the uth symbol which is equally likely to 
be ±1. The transmitted information bits are first encoded by 
a convolutional encoder. The rate of the encoder is R = 1/2. 
We have the code symbols c k , k € [0, 2U — 1]. Let the vector 
Cm = [ Co, ci, c m ] denote the first m(m < 2U — 1) code 
bits. We assume the modulation type is binary pulse amplitude 
modulation (PAM). The transmitted signal can be written as 



2(7-1 



s(t) = 2J c k p tr (t - kT s ) 



(1) 



fc=0 



where p a is the shaping pulse, and T s is symbol duration. 

The wireless channel can be modeled as a linear time- 
invariant (LTI) system with a finite impulse response h(t) = 
Yld=o oti&(t~Ti). In the case of wireless time-varying channel, 
we assume that within the coherent interval the channel can 
be modeled as time-invariant. The multipath delay spread of 
the indoor channel is T max . The ISI appears if T max > T s . 

The received signal can be written as r(t) = s(t) * 
h(t) + n(t), where * denotes convolution, n (t) is AWGN 
with double-sided power spectral density No/2. The system 
function of the ideal LPF is 



1/VWB \u\ < B 



others 

The variance of noise will be normalized to one after filtered 
by the LPF. The filtered r(t) can be expressed as 

2(7-1 

r(t) = Ck P^ ~ kT ^ + n ' W (2) 

fc=0 

where p K f(t) — p tI (t) * h(t) * h\ p (t), n'(t) = n(t) * hi v (t). 
hi p (t) is the impulse response of the LPF. 




Fig. 2. Joint decoding diagram 



Define sampling period T = \/(2B) = T s /N, which means 
for every pulse in the duration T s , N samples are generated. 
The received signal r(t) is sampled and quantized to one bit 
resolution. Let r m> j denote the ith sampling point in the mth 
symbol duration. We have 



+1 r(kT s +iT)>0 
-1 r(kT s +iT)<0 



(3) 



We can get the probability 

rn 

P(r m ,i = +l|cm) = Q(J2 c k p Kf {mT s - kT s + iTj) (4) 

k=0 

and 

m 

P{r m ,i = -l|c m ) = 1 - Q(y^ c k p Kf (mT s - kT s + iT)) (5) 

k=0 

where Q(x) = 1/\/2tt f °° exp(— t 2 /2)dt is the Gaussian Q 
function. 

The digital receiver can get r m ^ from the monobit ADC. 
The main work of the digital processing unit is to take the 
sampling point r TO ,j and estimate dk directly. 

III. Viterbi Decoding and DFE 

A. Joint Viterbi Decoding and DFE 

Joint Viterbi decoding and DFE is a method that estimating 
the original transmitted information bits dk- In order to make 
the BER as small as possible, we make the decision d k 
by its maximum likelihood (ML) probability. Fig. [2] shows 
the processing of the received points. In 0, an iterative 
approach was proposed to estimate channel state information 
(CSI) by transmitting training symbols. In our research, for 
simplicity, we assume the full CSI is known in advance. For 
the convolutional error correction decoder, we use Viterbi 
algorithm. 

Since the Viterbi algorithm stores information for each state, 
the complexity of the decoder is proportional to the number 



of states in the trellis. For a convolutional encoder that has 
fi registers, the number of states for this finite state machine 
(FSM) is 2^. Define vector r m = [ r m ,o, r m,i, -r m ,N-i } 
and matrix R k = [ ro, r±, ...r2k-i }■ In this algorithm, we 
need to construct a reference ISI waveform. 

2k— 1 

Pwav(*) = X] C "iPrAt - mT s) 

712—0 

c m is the estimation of c m . Similarly, Wm,i = pwav 

(mT s + iT) 

denotes the full resolution value of the ith 
point in the mth symbol duration. The vector 
w m = [ w mfi , w m ,u ».«;m,JV-i ] stands for N points 
in the mth duration. The matrix Wf. — [ wq, Wi, ... w 2 k-i ] 
contains the sampling points from the constructed waveform 
in the duration of ~ 2kT s . We choose proper W k so that 

d = arg max P(R k \W k ) 

At each step, we update matrix Wk- At the end of fcth 
decoding, cq, c%, ...C2k-i are decided. So for the (fc + l)th 
step, the ISI of previous symbols has been estimated. The 
current signal detection can be viewed as making a decision in 
AWGN channel. Thanks to the memorylessness of the AWGN 
channel, we have 



\ogP(R k \W k ) = logn^oVrfclu;*) 



(6) 



Every sampling point is independent with others. Therefore at 
the fcth step the log-likelihood probability is 



2k-l N—l 

\ogP{R k \W k ) = \ogp{r m 4w m ^) (7) 

Now we only consider the effect of noise. Then 

r m , 4 = sgn (w m s + n m<i ) 

where n m: i ~ N(0, 1). Define e m ^ = Q{w„ l: i), the probabil- 
ity of r m ^ can be expressed as follows: 

p(r m ,i = +l\w m .i) = 1 - e m ,i (8) 

p(r m ,i = -l\w mii ) = e m ,i (9) 
Combining <(8J an d ®, we get 

p(r m ,i\c m ,i) = 1/2 + r m>i (l/2 - e m)< ) (10) 

We substitute ( TTOb into (0, the log-likelihood probability is 
given by: 

2fc-l N—l 

log P(R k | W k ) = log(l/2+r m>i (l/2-e m>i )) (11) 
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Fig. 3. Decoding processing in the trellis 



B. Algorithm Description 

The joint Viterbi decoding algorithm attempts to find the 
maximum log-likelihood probability for each state. It gives 
the decoding result related to the current probability. As in 
the standard Viterbi algorithm, this algorithm is based on the 
trellis. The states and trellis are given by the convolutional en- 
coder. The algorithm needs to store the following information 
for each state. 

1) The surviving path leading to the state. 

2) The metric output of this path for the fcth step. That is, 
the Viterbi decoder's output for the "edge" from previous state 
to current state in the surviving path. 

3) The constructed ISI waveform corresponding to the 
surviving path. This waveform will be used by the next step. 

Initialization of the algorithm : 

Usually, we start the Viterbi decoding from the all-zero state 
So. The metric of so is set to be a large positive number. 
The algorithm works as follows : 

1) Construct a possible ISI waveform for each state. For 
coding rate R = 1/2, in the fcth (fc > 1) step, we get v^k-i 
and r%k-i from the received waveform in the duration of 2T S . 
We can construct 2 2 kinds of possible overlapped waveform 
related to {S 2 fc- 2 , ^fc-i} = {0, 0} , {0, 1} , {1, 0} or {1, 1}. 

2) For each arrived state, there are several leading paths. We 
compute log-likelihood probability for each path using (11). 
We also save the constructed waveform of the surviving path. 

3) Store the output of the current surviving path and its 
corresponding code bits. Make a decision of d k for each 
state. By choosing the surviving path, we make a decision of 
{c 2 fe-2,c 2 /c-i} for each state. Then we save {c 2 fc-2, c 2k -i} 
and d k . We also save the constructed waveform for each state. 

4) When the last step is finished, only one surviving path 
is remained. We need to trace back and get the final decoding 
result d. 

The joint decoding algorithm requires an amount of memory 
that is proportional to the number of states in the trellis. It also 
needs memory spaces to store the constructed waveform. The 
length is T max + T s . 
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Fig. 4. Path comparison 
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Fig. 5. 4 states 



C. Example 

Consider a typical trellis with the generator polynomial 
matrix g = [5 7] which is shown in Fig. [3] In the beginning, 
the metric of sq is set as M step o(so) = I, where I is a real 
positive number and / ^> 0. 

Fig. |3] shows the decoding processing in the trellis. For the 
first step, we receive Vq and r\, then our task is construct 
"possible" ISI waveforms for the arrived states sq and S2- 
There is only one coming path for each of them. We build 



and 



Plleplit) = -Pref(t) - Pref(t - T s ) 



P?tLl\(t)=Pref(t)+Pref(t-T s ). 



Pstepi(t) represents the constructed waveform for the path 
from so to S2 in the first step. Let M*° e pJ denotes the path 
metric from sq to S2 in the first step. The vectors Wq and Wi 
are sampled from p s s °^i(^) m tne duration of ~ 2T S and 
Mstepi i s computed using (fTTT l. We can get iu and Wi from 
Pltepi (*) an( i compute Mstepl m the same way. After the first 
step, we get 



M stepl (s ) = M step0 (s ) + M s J; s p l 



and 



M 



stepl 



(s 2 ) = M step0 (s ) 



We store the ISI waveforms, {c"o, ci}, do and metric for each 
state. The second step is similar to the first one, let's take the 
arrived state so as an example, we have 

PTtepKt) = Pfteflit) ~ Pref(t - 2T.) - Pref (t - 3T S ) 

From the third step, there is a competitive path for each state. 
In Fig.|U we can easily find that there are two paths leading to 
So- Under the assumption that {64, 65} = {—1, —1}, ^3 = —1 
for path 1, we build 

P?tepl(t) = PTtepUt) ~ Pref(t - 4T S ) - p ref (t - 5T S ), 

get and w 5 from p s s ° t e P ° 3 (t), and compute M^g. We 
do the same operation for the second path thus is 
achieved. The accumulated path metrics M step 2(so) + 
and M ste p2{si) + M^'g should be computed. We pick out 




Fig. 6. 8 states after expansion 



the larger one from the two accumulated path metrics, set it as 
the updated M ste p3(so) and denote its corresponding path as 
the surviving path. Finally, we store the constructed ISI wave, 
{64,65}, dz and M step3 (so) for the surviving path. Similarly, 
for the rest arrived states, we get surviving paths and store the 
information for each of them. When the last step is finished, 
we pick out the final surviving path which has the largest state 
metric and complete traceback. Finally, the stored dk in the 
final surviving path are the decoding result. As is shown in 
Fig. [3] if the final surviving path is so — > S2 — > si — > so, the 
result is d = {1, 0, 0}. 



D. State Expansion 

For a determined trellis, the number of states is given. 
The BER performance may be poor for lack of states. If we 
keep more states in the algorithm, we keep more interference 
situations that may arise. It is desirable to expand states. 
However, considering the algorithm complexity, we could not 
expand too many states. The amount of states will change after 
expansion. Consequently, the trellis also changes. 

Fig. [5] and Fig. [6] show an example of the original trellis and 
its expanded trellis respectively. The convolution code rate is 
1/2, generator polynomial matrix g = [5 7]. 

IV. Simulation Results 

The joint Viterbi decoding and DFE algorithm based on 
monobit digital receiver are applied to UWB signal processing 
in this section. Here, the second derivative Gaussian pulse [9| 
is used, which can be expressed as follows: 



Pa (t) = (1 - 4vr(Vr) 2 )exp(-2^(i/T) 2 ) 



(12) 



The constant r determines the pulse duration. 

The simulation conditions are as follows: 

We use standard CM1 multipath fading channel with 100 
realizations (8). The transmitted signal pulse is as (12) with 
r=0.22ns. The bandwidth of LPF is B = 5GHz. The sampling 
period is T=0.1ns. Each frame contains Nd — 2000 data. The 
transmitted symbol rate is 1GHz, that means T s =lns. We have 
the full CSI. The T max is about 100ns. Fading channel without 
ISI is also in use. For this type of channel model, we simply 
let T s =100ns. In this simulation system, channel coding is 
required. We use a typical convolutional code with the rate 
R=l/2. g = [171 133] is the generator polynomial matrix. 
The SNR is defined as E b /N = £ili P r 2 ef (iT). 

We explain the abbreviation in the simulation results. MB 
represents the monobit sampling. S represents for soft-decision 
Viterbi decoding. J represents joint decoding and DFE al- 
gorithm. The word CAS indicates that the separate DFE 
demodulation and Viterbi decoding method is used. The result 
for fading without ISI channel is denoted by NOISI. 

Fig. [7] shows the BER performance of the joint decoding 
algorithm and DFE- Viterbi concatenated algorithm in standard 
CM1 channel. 

MB-J-128S: Monobit digital receiver using the joint decod- 
ing algorithm that contains 128 states. 

MB-CAS: The DFE processing expands its state number to 
64 and Viterbi algorithm has 64 states. 

FR-S-NOISI: Full resolution detection and soft-decision 
decoding in CM1 channel without ISI. 

In Fig. H, we notice that the MB-J-128S has 3dB SNR loss 
compared to FR-S-NOISI. When SNR is larger than lOdB, We 
can observe about ldB gap between the two monobit curves. 
Fig. [8] shows the effect of state number. It can be seen that 
the BER decreases as the amount of state increases. The curve 
MB-S-NOISI is our final destination. As the number of states 
becomes larger, the gap between BER performance in the ISI 
channel and that in the channel without ISI becomes smaller. 
However, when the number of states increases to some extend, 




Fig. 7. Comparison of performance of joint decoding and separate decoding 
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Fig. 8. BER performance of joint decoding algorithm with different amount 
of states in fasing isi channel 



there is a BER floor. In the simulation result, state number 
changes from 128 to 256, the performance does not improve 
too much. By state expansion, we achieve about ldB gain. 
There is still about a ldB gap between the no ISI BER and 
joint decoding BER. 

V. Conclusion 

We presented the system model of the monobit digital 
receivers in ISI channel. The joint Viterbi decoding and DFE 
algorithm for monobit digital receivers was proposed. We gave 
the algorithm description and derived the log-likelihood proba- 
bility of the received signal. State expansion could enhance the 
BER performance of the algorithm. The BER of the separate 
decoding method and joint decoding approach was simulated. 
The simulation shows that the algorithm we proposed has ldB 
SNR gain. However, compared to the BER performance of no 



ISI channel, the algorithm incurs about ldB loss. There is a 
3dB gap between the joint decoding algorithm and FR-S-SEP. 
This method can be applied to UWB communication systems. 
It is also an efficient way dealing with ISI for monobit digital 
receivers. Future research topics include suppressing ISI with 
estimated CSI by training symbols, increasing the code rate 
R, improving the BER performance. 
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